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Abstract. We have used the radiative transfer code HDUST to analyze and interpret 
the long-term photometric behavior of the Be star w CMa, considering four complete 
cycles of disk formation and dissipation. This is the first time in which a full lightcurve 
of a Be star was investigated and modeled including both disk build-up and dissipation 
phases. Based on the quite good fit of the observed data we were able to derive the 
history of stellar mass decretion rates (including long- and short-term changes) during 
the disk formation and dissipation phases in all four cycles. 


1. Introduction 

to (28) CMa (HD 56139, HR2749; B3Ve) is one of the most observed southern Be 
stars. Long-term photometric monitoring exhibits a quasi-regular cyclic variation with 
an amplitude of about O.^S in the E-Band. This allows to study evolution of the disk in 
detail. In each cycle a new disk is formed during 2-3 years, and then dissipated in 5-6 
years. 

Carciofi et al. (2012) used to CMa to model, for the first time, the light curve of a 
Be star based on the viscous decretion disk (VDD) model (the theoretical description of 
the model is given by Haubois et al. 2012). The model of the dissipation curve allowed 
the authors to determine that the viscosity parameter of Shakura-Sunyaev (Shakura & 
Sunyaev, 1973) is a = 1.0 ± 0.2. Such value for a suggests that viscosity may be 
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produced by a disk instability whose growth is limited by shock dissipation (Carciofi et 
al., 2012). 
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Figure 1. Top: Observed light curve (grey triangles) of to CMa compared with 
our best-ht theoretical model (blue line). Bottom: Residuals of the ht. 
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Figure 2. The mass decretion rate history of to CMa: The red lines show the 
outbursts and the blue lines show the quiescence. 

The work of Carciofi et al. (2012) represented just a first step towards understating 
the physical conditions of to CMa’s disk. As such, the study had several limitations, 
among which we cite the fact that a was determined only during the phase of disk 
dissipation and by fitting only the V-band light curve. Lifting these limitation is the 
main motivation for this work. 


2. Results 

We used the time-dependent hydrodynamics code SINGLEBE (Okazaki et al. 2002; 
Okazaki 2007) to model the disk surface density evolution given a prescription for the 
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Figure 3. Consistency of model based on a =1.0 with the observed data even in 
the small events: This small outburst (phase I in Fig. 4) happened in the middle of 
main outburst phase of fourth cycle, after a small quiescence. Theoretical lightcurve 
is shown by blue solid line and the observed data is shown by the grey triangles. 


Epoch 

2009 2010 2011 2012 2013 2014 



55000 55500 56000 56500 

t (MJD) 


Figure 4. Theoretical lightcurve of co CMa for different viscosity parameters: a 
=0.1 (red solid line), 1.0 (blue solid line, 2.0 (green solid line: The lightcurve for a 
= 1.0 has broad consistency with the observed data, but does not fit well on the last 
rapid decline (Phase //). 


disk feeding history. To do that, several outbursts and quiescence phases with various 
time-scales, few weeks to few years, were considered. The surface density converted 
to volume density then is used as input for the three-dimensional Monte Carlo radiative 
transfer code, HDUST (Carciofi & Bjorkman 2006) to obtain the emergent spectral 
energy distribution (and other observables of interest), which is used to calculate the 
V-band excess, AT , of the disk as a function of time. 

For the results of Fig. 1 we adopted as stellar parameters the ones used by Carciofi 
et al. (2012) and let the mass decretion rate vary in amplitude with time, in a fashion 
described in Fig. 2. 
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One of the features of co CMa’s lightcurve is an obvious decline in the brightness 
of the system in the successive dissipation phases (green line in Fig. 1). We found that 
one possibility to interpret this behavior is to consider that at quiescence the star shifts 
to a lower mass decretion rate instead of entering a state of absolute quiescence. 

Our model shows the value of a should not be far from 1.0 at least for the most 
part. In Fig. 3 we show the first model of a Be disk build-up (for phase I in Fig. 4), and 
this model is broadly consistent with a alpha of unity. 

The dissipation phase marked as II in Fig. 4 presented a challenge for the model, as 
it is seemingly not consistent with and alpha of unity. Instead, a much lower alpha=0.1 
seems to better reproduce the data. This fit is still under investigation, since the slower 
dissipation in phase II may simply be the result of small-scale outbursts during the disk 
dissipation phase. 

According to the obtained results, we found that the mass decretion rate of co CMa 
displays a complex history including both short and long time-scales, also very different 
range of variation (Fig. 2). 


3. Discussion and Perspectives 

We present a complete model of the disk evolution of co CMa for the past forty years. 
For the first time, both outburst and quiescence phases were modeled with a single 
model. 

The overall fit is quite good, showing strong evidence of the VDD scenario as the 
mechanism controling the disk formation and evolution. 

A very puzzling feature of co CMa is the slow secular fading shown in the past 
forty years, in which each quiescent phase is less bright than the previous one. Here, 
we present one possible explanation for this phenomenon, by assuming that the mass 
decretion rate did not go to zero at quiescence, but instead assumed a smaller value 
(Fig. 2). The validity of this assumption needs to be ascertained by modeling other 
observables in addition to F-band magnitude. This will be the subject of future work. 
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